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A solution to the cosmological constant problem 
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We argue that, when coupled to Einstein's theory of gravity, the Yukawa theory 
may solve the cosmological constant problem in the following sense: The radiative 
corrections of fermions generate an effective potential for the scalar field, such that 
the effective cosmological term Acfr is dynamically driven to zero. Thence, for any 
initial positive cosmological constant Aq, Aeff = is an attractor of the semiclassical 
Einstein theory coupled to fermionic and scalar matter fields. When the initial 
cosmological term is negative, Acfj = Aq does not change. Next we argue that 
the dark energy of the Universe may be explained by a GUT scale fermion with a 
mass, m ~ 4.3 X 10-*^^ GeV {-sfko / IQi^^ GeV)-*^/^. Finally, we comment on how the 
infiationary paradigm, BEH mechanism and phase transitions in the early Universe 
get modified in the light of our findings. 



1. INTRODUCTION 

>< 

\ The cosmological constant problem can be stated follows Vacuum fluctuations of any particle 
species contribute to the energy density of the vacuum as, 



where p is the momentum of a particle, m its mass, and the positive (negative) sign stands for bosons 
(fermions). (In this article we use the natural units in which the speed of light and the reduced 
Planck constant are set to unity, c = 1 and % = 1.) po is a badly divergent quantity. It is often 
argued that a natural ultraviolet cutoff is the Planck scale, mp = GJj^^^ ~ 1.22 X lO^^ GeV, such 
that, 



mi 



of equivalently, 



Po = ±-P^ ^ ±(3.4 X 10^« GeV)^ , (2) 



Ao = ±-% - ±2.4 X 10=^^ GeV^ , (3) 
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where Mp = (8710^)'^/^ ~ 2.4 x 10^^ GeV is the reduced Planck mass. Above the Planck scale, 
an as of vet unknown quantum gravity theory is believed to lead to finite dynamics. The observed 
value 3,0] is at least 120 orders of magnitude smaller, 

Aobs < SH^n^E ^ 5 X 10-«^ GeV' , (4) 

where we used the current value for the Hubble parameter, Hq 1.5 x 10^^' GeV, and Qbe — 0.74 
is the dark energy contribution. 

The discrepancy between and (jH)) is known as the cosmological constant problem. There is 
a milder formulation of the problem, according to which supersymmetry takes care of the balance 
between the positive and negative contributions in Q above the scale at which supersymmetry is 
restored. In this case, po ~ i'^susy where msusy is the energy scale above which supersymmetry is 
restored. This works only for global supersymmetry however, since in local supersymmetry (which 
also includes gravity) this balance does not necessarily hold. 

Zeldovich has proposed |4] that the discrepancy between the observed (jH) and the theoretically 
more natural value Q could be explained by a mechanism, where quantum matter fiuctuations (j21) 
would cancel an initial geometric cosmological constant of Einstein's theory. Yet neither he, nor 
anyone else, has been able to provide a satisfactory compensatory mechanism. In this article we 
present such a dynamical mechanism, which explains how an initially large and positive Aq can be 
driven to zero by the gravitational backreaction induced by fermionic quantum fiuctuations, such 
that today the effective cosmological term has a small, but nonvanishing, value. 

The cosmological constant has a long history 0] since it was in 1917 introduced by Einstein 
to his general relativistic theory of gravitation. Around the same time Nernst observed that 
quantum matter (light) fiuctuations can contribute to the vacuum energy. In 1968 the vacuum 

n 

energy was experimentally observed j2| based on the Casimir effect Here we just briefly mention 
several attempts to solve the cosmological constant problem. It is well known that (tree level) scalar 
potentials {e.g. quintessence models j2|) cannot do the job because of the Weinberg theorem [l|. This 
theorem states that one can always add a constant potential energy without dynamical consequences 
to any scalar fleld potential which naively drives the fleld towards the value, where its potential 
energy vanishes. Another possibility is "shadow" matter, which contributes with an opposite sign to 
the energy density as ordinary matter, and couples only gravitationally. This requires introduction 
of a lot of new and yet unobserved particle species, and moreover has problems when perturbative 
gravity corrections are taken into account lul. The general trend in literature is either to resort to 
the anthropic principle |ll], or to postulate new symmetries, an example being the symmetry which 
relates positive and (as of yet unobserved) negative energy states |lol. 111! 11^. 
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When written in the approximation of a homogeneous and isotropic space-time with the metric 
tensor, 

5f^i. = diag(-l,a^a^a^) , (5) 
the semiclassical Einstein equations can be written in the following FLRW form, 



H = {-) = P + ^ 2 6 

a-' 6 6 a'^ 



d AuGn , An 

- = ^P + 3P+^, 7 

a 6 6 

where a denotes the scale factor, a = da/dt, d = d'^a/dt'^, H is the Hubble parameter, Gn the Newton 
constant, Aq the cosmological term, p the energy density, V the pressure, and k is the curvature of 
spatial sections of space-time. Here we take for simplicity = 0, which is consistent with current 
observations In the derivation of Eqs. (jHHZI) one assumes that the matter stress-energy tensor 
has an ideal fluid form, and Eqs. (jUHZI) are written in the fluid rest frame. 

We shall now calculate the stress-energy contribution to the Einstein equations (jBHZI) in the Yukawa 
theory. In doing so we shall resort to certain approximations, which are needed for analytical calcu- 
lations. 

The dynamics of the matter fields are described by the following model Yukawa theory, 

5'yu = j d'^X^f^CYn , (8) 

with the Lagrangian, 



where y = e|^7^(9^ — F^) denotes the covariant derivative acting on spinors, is the tetrad field, 
is the spin(or) connection, mo is the bare fermion mass, if) and (p denote the fermionic and scalar 
fields, respectively, Vq = (l/2)/iQ(y9^ + AoV5^/4! is the (tree level) scalar potential, Pq and Aq are the 
bare scalar mass and quartic self-coupling, respectively, IZ is the Ricci curvature scalar, g = det{g^^), 
and /o is the Yukawa coupling. 

Note that fermions appear quadratically in the action, so in principle they can be integrated out. 
In practice though, this can be done only in a handful of gravitational backgrounds. Since we are 
interested in solving the problem in a background given by a (positive or negative) cosmological 
term plus matter contribution, we make the first of our crude approximations and assume that the 
gravitational background is well approximated by (anti-)de Sitter space-time. Below we comment on 
how good this approximation is. The problem of integrating fermions in de Sitter background has 
recently been solved by Miao and Woodard 13], whereby they neglected scalar field fluctuations. This 
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approximation is accurate to leading log, where the log refers to ln(a), and a is the scale factor of the 



Universe (see Eq. fll4|) below). The leading log approximation is accurate in de Sitter background 15 1. 
Certain aspects of the work in Ref. [l^ are inspired by an earlier work of Candelas and Raine 16 1. 



Various quantum aspects of the Yukawa theory in de Sitter background have former 



in Refs. 



Ulll^ll^; and a new fermion mass generation mechanism has been found 12 1. 



1. 



Deen studied 



2. DE SITTER AND ANTI-DE SITTER SPACES 

de Sitter space is perhaps best viewed as a 4-dimensional hyperboloid embedded into the five- 
dimensional Minkowski space-time with the line element, 

dsl = -dXl + dXl + dXl + dXl + dXl . (10) 

The embedded hyperboloid of de Sitter space is shown in figure [H and it is determined by 

-Xl + Xl + Xl + Xl + Xl = ^^, (11) 

where H denotes the Hubble parameter. The symmetry of de Sitter space, S'0( 1,4), is manifest in 
this embedding. One can define de Sitter invariant distance functions, 

^ 1 

Z{X-X')=H^ J2 VabXaX'^ = 1~-Y{X;X'), r/^s = diag(-l, 1, 1, 1, 1) . (12) 

A,B=0 ^ 

We shall use the following flat four-dimensional coordinates (which cover 1 /2 of de Sitter manifold) , 

1 H 
Xq = — sinh(iJt) + — e^*||ff , (-oo < t < oo) 
H 2 

Xi = e^^Xi (-CX) < Xi < oo) , (« = 1, 2, 3) 

X4 = 4 cosh{Ht) - -e^*||f ||2 , (13) 
H 2 

in which the metric tensor reduces to the form (0), with the scale factor, a = e^*. When written in 
terms of conformal time 77, defined as adrj = dt, the metric tensor acquires the conformal form, 

9i,u = a?r]^^, a = --^ (^ < 0) , 77^^ = diag(-l, 1, 1, 1) . (14) 

The invariant distance functions, Z{X;X') = z{x;x') and Y{X]X') = y{x;x'), reduce in these 
coordinates to the simple form, 

z{x; x') = I — -y{x; x') , y{x; x') = aa'H'^Ax^ (15) 

with a = a{ri) = —l/{Hri), a' = a{rj') = —\/{Hrj'), and 

Ax^(x; x) = — T]'\ — ieY + ||x — x'H^ , (16) 
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FIG. 1: The embedding of de Sitter space into a five dimensional flat space-time. The vertical line corre- 
sponds to the time coordinate, Xq = T, and the radial coordinate, R = {Xf + Xl + + X^y^'^ = 1/H. 
At each point {T,R) there is a unit 3-sphere, (which has in this figure collapsed to S^). 

where (for later use) we introduced an infinitesimal e > 0. In these coordinates the curvature of 
spatial sections vanishes, and thus they are also known as flat (Euclidean) coordinates, in which de 
Sitter space appears as uniformly expanding. 

For the integration of fermions we also need the fermion propagator. 



iS{x; x') = {x\ 



\x') , (m = mo + fo(p) 



(17) 



This propagator can be solved in de Sitter background and in the approximation of a nearly constan ; 
ipid^Lp and higher order derivatives are neglected), and the solution can be written in the form 



iS{x; x') = a(i^+ {aa) 2 ^ 'iS±{y) 



(18) 



where = a 2 i^'^df^a'^^ and iS-i- = iS±{y) are the following de Sitter invariant scalar functions 
(which contain no spinor structure), 

D ■m\j^ I D 

~r 1.— \ .11 n-Y^ I I /no II n, . 

(19) 



where m = m{ip) is defined in Eq. (jl7p . D denotes the number of space-time dimensions, 2-^1 is the 
hypergeometric function, and 7^^ are the Dirac 7-matrices. With the e pole prescription as in (|15flTH|) . 
one can show that the propagator (fTSHT^ indeed solves the Dirac equation in de Sitter space, 

^/^{^-m)iS{x;x)=i5'^{x-x'), (20) 

where 5^{x — x') denotes the D-dimensional Dirac (5-distribution and (p = const. 
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Similarly to de Sitter space, anti-de Sitter space can be thought of as the hyperboloid, 

- -Xl + Xl + Xl + Xl = -R^^-^^, (21) 

embedded into the five- dimensional fiat space-time with the line element, 

dsl = -dXl - dXl + dXl + dXl + dXl , (22) 

which possesses an SO(2,3) symmetry. The graphical representation of anti-de Sitter space is the 
same as in figure ^ except that each point on the vertical axis corresponds to a circle, and the points 
of constant radii, R = (X| + X^ -|-X|)^/^, correspond to two-dimensional spheres 5^. The definition 
of anti-de Sitter invariant distance functions is as in de Sitter space (|T^ . with the five dimensional 
metric given by, tiab = diag(— 1, —1, 1, 1, 1). 

The following coordinates, which cover a 1/2 of anti-de Sitter space, are analogous to the fiat de 
Sitter coordinates (fT^ . 

Xo = ^ + ^(^R^ + xl + xl-r]'^) (-cx)< r/ < oo, n > 0) 
Xi = RuT] 

Xi = RuXi {i = 2, 3) (— oo < Xi < oo) 

X, = l--^{R^-xl-xl + v'), (23) 
With the additional replacement, u = 1/xi, these coordinates yield the conformal metric tensor, 

9f,u = a'^Vi^u , a = {xi>0), ^ ^ ^'^^^ 

One should not confuse H in anti-de Sitter space (where it denotes the inverse radius of curvature of 
the space) with the Hubble parameter H in de Sitter space. Note that unlike de Sitter space, which 
represents an expanding space time in conformal fiat coordinates, the conformal fiat coordinates of 
anti-de Sitter space correspond to a static slicing with the space-time curved in one spatial 
direction, which we choose to be Xi. One can easily calculate the anti-de Sitter invariant distance 
function, Z = —H'^J2a,bVabXaX'b = 1 + Y/2. Upon defining, y{x;x') = Y{X;X') and z{x;x') = 
Z{X;X'), we find in these coordinates, 

z{x; x') = 1 + -y{x; x') , y = aa'H^Ax'^ , (25) 

where here a = a{xi) = l/{Hxi), a' = a(xi') = l/{Hxi') {xi,Xi' > 0) and Ax^ is given in Eq. (fT^ . 

One can show that the anti-de Sitter fermion propagator in the presence of an approximately 
constant scalar field reduces to the form [2^, 

1 ±^7^ 

± 



iS{x;x') = a{if+rnj{aay^J2'^S±iy)-—f^^ (26) 
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where iS± = iS±{y) are the anti-de Sitter invariant scalar functions, 

^ ,<--../^i^^fcl|M±l) ,,J£.^^^,£^^Al^^, (27, 

^ (47r)f r(f) ^2 H' 2 H' 2' 4^ ' ^ ^ 

and m = mo + fo(p. Note that the two propagators (fTHHT^ and (jlHHlZI) are related by the analytic 
continuation, 

H-^iH. (28) 

The extra phase in (j27|) comes from the term in the propagator proportional to y^'^/"^. This term 
gives rise to the 5-function in the propagator equation ((201), which induces an extra phase due to, 
y = aa'H'^Ax'^ —y. The spinor structure of the two propagators is related such that the positive 
and negative energy projectors are replaced as, (1 ± 7°)/2 ^ (1 ± i'~f^)/2. The extra imaginary i in 
the projectors, (1 ± i7^)/2, is important, since it assures that the spinor structure drops out from 
the differential operator acting on the propagator. 



EFFECTIVE POTENTIALS 



When one uses the anti-de Sitter propagator (|26H27|) to calculate the effective action by an anal- 
ogous procedure as done in Refs. and one obtains the following (renormalised) effective 



potential, 
where, 



El 

87r2 



''m\2 , ^ ^ /m\4 r^/H 



(29) 



Kree(^) = \f^V + ' , (30) 



ip{z) =^{d/dz) ln(r(z)) is the digamma function, and A, /i and ^ are the renormalised parameters (in 
Ref. Q] the renormalisation has been chosen such that these parameters vanish), ^(3) ^ 1.202 is the 
Riemann zeta function, and 'Je — 0.577 is the Euler constant. This effective potential is induced by 
fermionic fluctuations in the presence of a nearly constant scalar field in anti-de Sitter space. Since 
the integrand has simple poles at x = 2, 3, 4, .., the lower limit of integration can be taken as Xq = 
for m/H < 2, and Xq = [m/H] + CAds when m/H > 2, where < CAds < 1 is a constant, and [■] 
denotes the integer part. For m/H > 2 the potential ()29|) contains an unspecified constant, which 
has no dynamical relevance (the dynamical relevance is associated with —dV/d(p). 

Just like the propagators, the effective potentials are related by the analytic continuation (j2Hl), 
such that the de Sitter effective potential reads [l^ . 



V{^) = \4ree(<^) + 



H 



'm\2 , ^ ^ /m\4 /■m/_ff 



-27s(-) +(C(3)-7ij)(-) -2 dx{x + x')(i;{l+ix) + ^{l-tx)) 



.(31) 
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The large field limit of this effective potential is, 

^^-^InQ (m»i/). (32) 

The potentials ()29fl!?T|l determine the dynamics of scalar fields in de Sitter and anti-de Sitter back- 
grounds. The dynamics in de Sitter space is given either by the corresponding Euler-Lagrange 
equations, 



dl^--d^^ + 'iHdt^ = -—, (33) 



1 dy_ 

nn ^ 

or by the Starobinsky stochastic theory |22| [1^, according to which 



where ^ is the white noise 



i^i^a^))^^Ht^e fff-:^ , (35) 

477"^ Ha\\x — X \\ 



generated in accelerating space-times as vacuum fluctuations get stretched beyond the Hubble radius, 
1/H. When the effects of quantum tunneling and stochastic fluctuations are neglected, the field 
generally evolves in the direction of decreasing potential, dV/dLp < 0. 

In anti-de Sitter space the dynamics is determined by the Euler-Lagrange equation, 

-^dlv - + 2H^^^ = (36) 

' a oxi dip 

where V is given in ()29|1 . Note that the "damping" (third) term in anti-de Sitter space generates 
inhomogeneities and thus has a different meaning from its de Sitter counterpart in Eq. ()33|) . Just like 
in de Sitter space, the Euler-Lagrange equation ()36p tells as that the field ip is driven in the direction 
of decreasing potential (along which —dV /dp is positive). 

In the large field limit (m = mo + /oV^ ^ H), the effective potential (1311) reduces to the form that 
deceptively looks like the Coleman- Weinberg effective potential Q, IigI. I21I , which exhibits the well 
known instability, according to which fermion fluctuations drive the scalar field to infinity, — > 00, 
whereby V ^ —00. This instability has been declared a problem. As we shall argue, when interpreted 
correctly, this runaway feature is the cure. 



4. STRESS-ENERGY TENSOR 



We have so far said nothing about the back-reaction of matter fields onto gravity, as indicated in 
Eqs. (JOHZI)- To address this question, we need to calculate the stress-energy tensor of the Yukawa 
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FIG. 2: Naive potential energies, v = vq + (167r^/i?^)V^, as a function of m/H = {niQ + fo^) /H for de Sitter 
(left plot) and anti-de Sitter background space (right plot). Note that, as (p increases, the potential energy 
decreases (increases) in de Sitter (anti-de Sitter) space. For definiteness, we chose ^ = 0, /x = 0, A = and 
u(0) = 100 (—100) in de Sitter (anti-de Sitter) space. 



theory (j^I), which is obtained by the variation of the action (jHJ with respect to the metric tensor (or 
tetrad). The relevant part of the stress-energy tensor is (here we suppress terms, which involve space- 
time derivatives acting on scalar fields, and whose contribution tend to be suppressed in accelerating 
space-times), 

Tf,u = ^\im^^^,(i'y^^V^)iS{x; x') - iS{x; x') V(^i7i.)) - ^oG^^ip'^ + ^/io0^ + ^0"^ + C((9v9)^) , (37) 

where 7(^Vy) denotes symmetrisation with respect to and u, and G^,^ = 71^^, — ^g^ulZ denotes 
the Einstein curvature tensor. Upon evaluating T^j, in anti-de Sitter background and performing 
dimensional regularisation and renormalisation, one obtains |2fl] (m = mo + /oV'), 



T. 



-Vs{v)g,. + o{{dvf) 



Vs{'p) = Vtree 



(3^ 



167r2 



(27.-l)y +(- + 2C(3)-27.)(-) +(-) (l-(-) + 



The de Sitter potential energy is, as expected, obtained by the analytic continuation p8|l of Eq. 



Vs{ip) = Vtree 



(39) 



167r2 



(2^--^)y -(2+2C(3)-27.)y +y (l+(^ 



m\ / m 



which agrees with the result obtained in . We plot these two potential energies in figure |21 Note 
that the de Sitter potential energy shows instability and decays for large values of the argument. 
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while the anti-de Sitter potential is stable. As the argument approaches, m/H 2, anti-de Sitter 
potential grows and diverges as, Vg +00, such that the two curves cross. For each point on the 
de Sitter curve, there is an anti-de Sitter curve which crosses it. If they were both correct, the 
field dynamics would drive the field both towards larger values (along the de Sitter curve) as well 
as towards smaller values (along the anti-de Sitter curve). Obviously, both cannot simultaneously 
correspond to the correct dynamics. 



20 

-20 
V -40 
-60 
-80 
-100 
-120 



AdS 
VdS 
AdS 







0.5 



1.5 

m/H 



2.5 



FIG. 3: Improved, but still naive, potential energies, Vg = vq + (167r^/ff^)V's, as a function of m/H = 
(mo + foip)/H for de Sitter and anti-de Sitter space-times. Two (anti-)de Sitter curves are shown. The 
initial potential energy for the de Sitter and anti-de Sitter case is chosen such that they match at Ves = 0. 
As in figure 121 we chose, ^ = 0, ^ = 0, A = 0. 



Peflf ~ V, 



eflf 



(40) 



Here Mp = {SttGn) ^ — 2.4 x 10^^ GeV denotes the reduced Planck mass. This leads to an improved 



This inconsistency is resolved by noting that de Sitter potentials are valid for all p^s > 0, while 
anti-de Sitter potentials are valid for pefr < 0, pefr = {H = 0) being the limiting case, where 

Vs + M^Ao , for Vcs > and Aq > 
Vs + M|Ao , for Kfr < and Aq < 

-1 ^ o /I XX 1 r)l8 

dynamics, two examples of which are shown in figure El Note that the anti-de Sitter potential energy 
has simple poles at m = nH, where n = 2,3,4, .. are the integers greater than or equal to 2, such 
that, for example, the second anti-de Sitter curve in figure El runs off to minus infinity at m = 2H. 
The first reaction of a learned reader may be: if this picture were true, that would be disastrous, 
since that would imply that Minkowski space with Aes = would be unstable, and decay into anti-de 
Sitter space, from where it would run either to some large and negative value of A^s (when m/H < 2 
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at Vcs = 0) or even worse to A^s — oo (when m/H > 2 at Vcs = 0). This results in a very 
uncomfortable observation that an unfortunate field fluctuation or a tunneling event might take us 
to anti-de Sitter space, resulting in dire consequences. There is one way we can be saved: if the field 
evolution slows down dramatically as one approaches Agg = from above, we would never reach 
anti-de Sitter. As we shall argue, this is indeed what happens. Even better, we shall see that the 
naive interpretation in figure |21 is incorrect. In fact, the de Sitter and anti-de Sitter curves do not 
continuously connect at Acfr = 0, as suggested by figure |2l and the curves with m > 2H are in fact 
unphysical. Moreover, Agff = is the stable end-point of the evolution for all Aq > 0. We shall also 
see that the probability to tunnel to anti-de Sitter space is vanishingly small. 



5. DYNAMICAL RELAXATION OF THE COSMOLOGICAL TERM 



The true dynamics in de Sitter space are obtained by taking account of the matter back-reaction 
via the Einstein equation When kinetic energies are neglected with respect to the potential 
energy Ves in (j^ . which is justified in de Sitter and quasi-de Sitter space-times (slow-roll regime of 
primordial inflation), the Friedmann equation (jHl) reduces to. 



3Mi 



1 / , Vs\ 



(41) 



ip 3 - ivip- 

with the scalar field dynamics given by Eq. or by Eqs. 

Taking account of the matter backreaction means simply that Eq. (jlT|l should be understood as 
H = H{(p). The true potential energy V^fj = Kfr(v^) (or equivalently the effective cosmological term, 
Aeff = Ves/Mp) is then given by the self-consistent solution of equations (|iT| and 

To get analytical insight into the field dynamics, observe that in the limit when m/H = (mo + 
fo(p)/H ^ 1, the potential energy (j39|l reduces to the form. 



Vs{ip) ~ Vtree " 



m\ 4 
H 



ln(-)-- + 7.-C(3) 



/m\2 



1 5 



+ 



11 

120 



(42) 



where Vtree is give n in Eq. (pUjl . At a first sight the potential seems to be of the Coleman- 



Weinberg type 



21j , whereby H is interpreted as a renormalisation scale. Moreover, for a sufficiently 



large expectation value of the field, the potential seems to be unstable and naively drives the field 
into anti-de Sitter space. Yet when Eq. (PT| is taken account of, the picture completely changes. To 
see this we insert (jl^ into (jUj) to get. 



Ao 
3 



m 



37rmp 



(m > H) 



(43) 
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where mp = Gj^^"^ ^ 1.22 x 10^^ GeV and we kept only the leading order term in (j^^ . If we want 
to solve this for H or Apff as a function of m, the following form of (j43p is suitable for iterative 
procedure, 

ff = mexp( ^ + ^— . 44 

When ^ m^/(37rm|) ~ (Ao/3)/ ln(m//f) (which is a stronger condition than m ^ if), already 
the leading order iteration becomes a good approximation, 

ff.,„exp(--£^) K«3^). (45) 



A careful look at (j45j) reveals something quite bizarre. The potential energy corresponding to (j45p is 
simply. 



Kes = Sm^M^xp ( - ^^^) , (46) 



27rmpAo^ 

which is of a completely different form than the perturbative-looking potential (021) • Moreover, the 
potential is nonperturbative with respect to both gravitational and Yukawa coupling constants. 
Furthermore, one can show that the full solution of Eq. ()43|) is also nonperturbative with respect 
to Aq. Recall that, after integrating the fermions in a nontrivial background (de Sitter) space-time, 
we obtained a seemingly perturbative effective potential. Yet, when the matter backreaction onto 
gravity is taken account of, the resulting effective potential becomes nonperturbative both in the 
Yukawa coupling, as well as in the gravitational couphng constants (the latter is quite common in 
classical gravity). 

The effective cosmological term as inferred from Eq. is, 

2 / 27rm|AoN 
Aeflf = 3m exp — , (47) 

which approaches zero if m approaches zero. Recall that m = mo + /oV^, such that this holds 
independently of the initial fermion mass mo. What remains to be shown is that m — at late 
times. It suffices to show that the potential energy pUj) (at this level of approximation, V = Vg, cf. 
Eq. ()32|)) drives the system towards m = at late times, where A^s = 0. Indeed, since within our 
approximations, —V = —V^ oc —m' = — /o < {V = dV/d(p), the force pulls the field towards the 
origin at m = 0, where both the potential energy Vg and = Acfr/3 exhibit essential singularity. 
Visually this can be seen in figure EJ in which we plot the effective cosmological term (j47p as a function 
of m. The arrows in figure |3] indicate the direction along which the field evolves. To conclude, we 
found that the cosmological constant problem is solved in the Yukawa theory plus Einstein's gravity 
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FIG. 4: The effective cosmological term Agg = 3H^ as a function of m = mo + fo^- The left pannel 
shows log^oC^eff), the right pannel AcS- As m decreases, Aeg approaches zero faster than exponentially. 
The dashed (blue) and solid (red) curves correspond to the 0th and 1st order iterative solutions of Eq. (|44() . 
respectively. From the left pannel we see that the effective cosmological term, A^s = 5 x 10~^^ GeV^ 
in Eq. Q, corresponding to the dark energy density of the Universe today, implies for the fermion mass, 
m = 4.3 X 10^^ GeV (where we chose V^o = 10^^ GeV). 

in the following sense: independently of the initial value of the cosmological term Aq, as long as it is 
positive (and smaller than the Planck scale), it is dynamically driven to zero. 



Aeff 0, 



(48) 



Hence Acs = is an attractor for any positive initial cosmological constant Aq > in a Yukawa 
theory coupled to gravity. This is true for any Yukawa theory whose tree level couplings fi, A and ^ 
are such that the scalar field does not get stuck forever in a false vacuum. 

Next we address the question of tunneling from de Sitter to anti-de Sitter space. The tunneling 
probability, Ptunneh can be estimated in terms of the Euclidean action Se as follows. 



tunnel 



-Se 



Se = J d'^XEy/gEMpAes 



(49) 



Working in closed coordinates (which cover whole de Sitter space), 
1 



^,2 



■(^dril + dx^ + sm\x)ide^ + sm\e)d(l)^ 



oo < TjE < oo, < X < tt) , (50) 



if 2 sinh (tie) 

we have ^/qe = sm'^{x) sin(6')/ smh.'^{r]E) such that the EucUdean action becomes, 

3^ 



^sinhV^)fsinh^(^) + - 



(51) 
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where t^; is the physical Euchdean time, which determines the instanton duration, and it is defined 
as, e^''^ — 1 = 1/ sinh^(ift^/2). Note that the action ()51|) diverges as t^; — > oo, implying that the 
tunneling is strictly speaking forbidden. For a Hubble time instanton, ~ Eq. ()51|) yields, 

Se-'§, (52) 

such that the probability for Hubble time instantons is suppressed as, Ptunnei ~ exp(— mp/iJ^), which 
is so tiny that it is indistinguishable from zero today. A typical time interval that anti-de Sitter 
instantons can exist (for which Se ~ 1) can be estimated from (j^T|) to be, Ie ~ (mpif)~^/^ ~ 2 x 
10^^ GeV^^ ~ 10^^^ sec. The question of the role of stochastic fiuctuations is under investigation |2o| . 

Having established that de Sitter space is stable against rolling and tunneling into anti-de Sitter 
space, we now comment on what happens if the Universe begins with a negative cosmological term 
Aq. The relevant potential energy curve is shown in figure El (the lower left [dashed green) curve). The 
field begins at some Aq < 0, and since the potential is concave (curved upwards), the field fiuctuates 
around m = 0, as given by Eq. (jHHI)- Just like in de Sitter space, the potential energy plotted in 
figure El is naive, and the full account of the potential energy (when field fiuctuations are suppressed) 
can be obtained by the self-consistent solution of Eq. ()38p and the corresponding Friedmann equation 
which can be obtained from Eq. © by analytic continuation. 



2^/1 da^2 _ _Ao _ Vsjif) , . 



where here H = 1/i? is the inverse radius of curvature of anti-de Sitter space, and we approximated 
the energy density by the corresponding anti-de Sitter potential energy (j^Hj) . p ~ K(. For a large 



curvature radius, R = H ^^y3/|Ao|, the potential can be expanded around the singular point 
m = 2H (which is approached whenever |Ao| ^ nip) to obtain (we set /x, A,^ = 0), 

"-■^^^-'-^(j^^^^'^-l) (^^«V3). (M) 
Just like in de Sitter space, the self-consistent solution in anti-de Sitter space exhibits two branches. 
The effective potential at the upper near-Minkowski branch (where H ^ 0) can be approximated by, 

Vs ~ -AoM^ - ^m^Ml , (55) 



such that 



Aeff = -3/j2 ^ -^m^, (56) 



which behaves like an inverted harmonic oscillator with a negative mass term (c/. Ref. |28l|). 

<ds ^ -^/o^p (^0 = 0, « -Ao/3) . (57) 
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In this analysis we have also assumed, m <^ (m||Ao|)^''^. These results can be interpreted as follows. 
If the Universe begins with a negative initial Aq and at the lower branch, then it is dynamically 
driven towards the stable point at, m = 0,AcfT = Aq, which is the true potential minimum. If, on 
the other hand, the Universe begins at the upper branch (j55p . it will exhibit an inverted harmonic 
oscillator instability with the mass (j37|) and quickly descent to the lower branch and eventually end 
up again at the stable point, A^g = Ag. In this sense there is no problem with stability of Minkowski 
space. 



6. DARK ENERGY 

Next we address the question whether the dark energy of the Universe , whose density is 

Qde = — = 0.74 ±0.04, (58) 

Per 

(the 3 year WMAP data [3| with a Hubble parameter prior results in, Qde = 0.76^g Qg) or equiva- 
lent ly 

AcfT,DE ^ 5 X 10"^^ GeV^ , (59) 

can be explained within our model, where, pcr = Si^pMl ~ 4 x 10^^^ (GeV)^ is the critical energy 
density today. To see that, we recast Eq. ()45|) in the following form suitable for iterations. 



4 vrmpAo /^^o 



from where we obtain, 



m 



m 



\n{m/H) ^ Snm 



{H' « ^) , (60) 



The choice Aq = (10^^ GeV)^ is the natural scale of primordial inflation |2^. Note that m oc 
Ag*^^ — H]^^ {Hi is the Hubble parameter during primordial inflation) changes rather slowly, as the 
scale of primordial inflation changes. Even though m in Eq. ()61|) is of the GUT scale, the curvature 
of the effective potential ()46|) today is of the order the Hubble parameter today, (PV/dip^ ~ Hq, 
which is tiny. The mass (1611) is our prediction for the heaviest fermion which acquires mass via the 
mechanism presented here. There may be heavier fermions in Nature, which couple to a (grand- 
unified) Higgs field [2^, e.g. to a field with a negative mass term and a positive quartic coupling, 
but they do not contribute to ^de- Indeed, the residual vacuum energy of all Higgs fields and 
other scalar and condensate fields is already contained in Aq, such that the dynamical compensation 
mechanism pres ented here works for the residual vacuum energies of any Brout-Englert-Higgs (BEH) 
mechanism |26| [2^ and phase transition, which takes place in the early Universe. What is special 
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about the scalar field responsible for our mechanism, is that it must not couple to gauge and other 

scalar fields which would result in an upward sloped effective potential. 

What remains to be shown is that, from the moment when the energy density becomes dominated 

by the residual fermionic energy density Vj-es in ()46|) . the Universe enters a slow roll regime. We thus 

need to estimate the rate of change of the Hubble parameter. From Eq. (j45|) we infer, 

H ^ 47rm|Ao m ^ 
H m 

In the slow roll regime we have, ip ~ —V'/{?)H). This and Eq. ()46p allow us to estimate m/m, such 
that (jU^ gives, 

H ^ 4vrm^Ag 

Slow roll regime requires, \H / H'^\ <^ 1, which means that 



« = TT^^JT- - X 10-%^^^) ^ (64) 



^mlAo ^1013 GeV' 

This is a rather small value, but still corresponds to a modest fine tunning. The electron Yukawa is 
sufficiently small, /e ~ 2 x 10~^, to marginally satisfy the bound (jMj) . 

From the rate of change of dark energy, pde — Vres, which is defined in Eq. (HHI), 

= -3/7(1 + wde) = -Svr/o , (65) 

we immediately infer. 



^Aq 1 mp 



/2 127r3/2 1 + ^DE 



ln( 



vrmpAo 



5/2 1 X 10^4 GeV 
1 + wde 



(66) 



When the slow roll condition (jMj) is met, wde gets close to —1, as it should. In evaluating 
we took, Ag^^ = 10^3 GeV. The r.h.s. changes by about 10% for each order of magnitude change 

1/2 

in the scale of primordial infiation, Ag . Since the current observational bounds on wde, —1.2 < 
('?i'DE)obs < -0.8 and (wDE)obs = — 1-06_„„„ do not resolve wde from —1, we cannot yet 

use (jUBj) to constrain the fundamental parameters, /g and Ag, of the theory. 

Next we estimate the number of e-foldings and the time before the Universe reaches the singular 
point at AcfT = (fiat Minkowski space). The number of e-foldings can be estimated as iV' = dV/dip), 

Since we do not know the value of /g/v^A^, we cannot yet calculate A^. Based on the slow-roll 
bound ()64|) we can however estimate the minimum number of e-foldings before the Universe hits the 
singularity at m = 0, 

2 /m.\ , , 
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The amount of time left before the Universe reaches the singular point can be estimated as, 

with uq = ('7rmpAo)/m'^ = In {m/H) ~ 132. Hence, the singular point Aegf = is reached only at an 
infinite future, which protects us from the singularity. 

Flat Minkowski space is thus the singular point of Yukawa theory, which is reached only at future 
infinity. If that point were ever reached, strictly speaking one would loose all information about the 
initial state, since all universes with an initial Aq > converge to that point, representing a new type 
of information loss. Figure 01 is thus incorrect in the following sense. If one ever reaches the point 
Acfj = 0, m = 0, any information about the initial Aq would be lost, such that it makes no sense 
to continuously connect de Sitter and anti-de Sitter curves, as indicated in figure El For a negative 
Aq, the evolution in anti-de Sitter space is trivial, in the sense that the effective cosmological term 
does not change, A^s = Aq. If for some reason the Universe starts at the upper near-Minkowski 
branch (jHU, Aefr ~ 0, then it will rapidly descent to Acfr = Aq < 0, as explained at the end of 
Sectional This is really an instability associated with anti-de Sitter space. Minkowski space is the 
space with Aq = 0, and it is thus fully stable. 

When viewed as a function of the fermion mass m, as m becomes smaller and smaller, the singular 
point Acff = is approached faster than exponentially, implying that at late times trajectories with 
different initial Aq rapidly bunch up, and as time goes on it becomes more difficult to resolve the 
original Aq of the Universe. Since our theory contains essentially three free parameters (Aq, /o and 
m), to resolve Aq and /o, in addition to ^de and wde, one has to measure dw^E/dt, which is planned 
to be measured in the near future. 



7. DISCUSSION 

A dynamical mechanism for relaxation of the cosmological term is presented in section Elby making 
use of the Yukawa theory coupled to Einstein's gravity. The mechanism is generic however and should 
apply to any theory, which in Minkowski space exhibits an effective scalar potential with instability 
or runaway behavior. When the matter backreaction to Einstein's gravity is taken into account, the 
scalar field is dynamically driven towards the attractor at Aeff = 0, for any initial positive cosmological 
term below the Planck scale, Aq ^ rrip, and for any reasonable tree level coupling parameters. If 
Aq < 0, the dynamics is trivial, and Acs = Aq. 

Related work on Yukawa theory does not offer an explanation for why Acfj ~ today. For example, 
in Ref. |16|] we read, "there is no particle creation in de Sitter space." In references which discuss 
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the Coleman- Weinberg-type of effective theories, it is universally claimed that the effective Yukawa 
theories pose a problem, because of the instability of the effective scalar potential induced by the 
radiative effects of fermions. 

An early attempt to attribute dark energy to the quantum fermionic fluctuations in expanding 
space-times dates a couple of years back . The attempt failed since at that time the authors were 
not aware of Ref. and did not properly calculate the contribution of fermionic fluctuations to the 
stress-energy tensor in de Sitter and other expanding backgrounds. 

It is worth noting that in the light of our dynamical relaxation mechanism for Agfj, the prob- 
lem of constructing quasi-stable local minima with a small but positive effective cosmological term, 
an example being the recent construction within the context of string theory 2^, loses its main 
motivation. 

The next question is whether it is reasonable to assume that the effective scalar potentials in de 
Sitter and anti-de Sitter spaces are dominated by the Yukawa coupling to fermions. It is well known 



that coupling to scalars contributes positively to the scalar potential energy |16l . 0| (in the latter 
reference, replace by A<y9^). It is less known that gauge fields contribute also positively to the 
effective potential js^, such that for the dynamical mechanism to be operative as advocated, it must 
be that the coupling of scalars to fermions dominates over the coupling to scalars and gauge fields 
at large scalar field expectation values at least for one scalar field. If the heaviest fermion of that 
kind has a mass, m ~ 4.3 x 10^^ GeV {y^Tio / 10^^ GeV)^/^, and the Yukawa coupling not larger than 
the electron Yukawa, then our mechanism can explain the dark energy of the Universe, which makes 
up about 74% of the energy density of the Universe, and which has a negative equation of state, 
wj)E e (-1.14,-0.93) 0- 

The history of the Universe gets revised within the Yukawa theory (jHHSI) as follows. The Universe 
begins with an initial positive cosmological term, Aq > 0, and with a scalar field vacuum expectation 
value close to zero. The latter can be achieved, for example, with a positive scalar mass term or quartic 
self-coupling. The initial cosmological term is not just of geometric origin, but it also comprises 
contributions of vacuum fluctuations of all matter fields. Even though the precise value of Aq is not 
important, if we want to realise primordial inflation within our model, then Aq ~ (10^^ GeV)^. If 
the initial Aq differs significantly from (10^'^ GeV)^, than one or a series of phase transitions in the 
pre-inflationary Universe can eventually generate a value suitable for primordial inflation. What is 
important is that the last slow roll regime of the theory corresponds to a scale, a/Aq ~ 10^'^ GeV, 
since then the amplitude of cosmological perturbations will correspond to the measured value. Since 
vacuum fluctuations of fermions in de Sitter background generate the effective potential ()31|) . the 
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field rolls driven by, ~dV/dip > 0. If the potential happens to be trapped at a small value of 
the field by e.g. a positive scalar mass term, the field eventually tunnels to a value above which 
—dV/dif > 0. One gets a slow roll inflation and enough of e-foldings when the condition, /o < 10~^, 
is met (note that this condition is milder than the slow roll condition today (jHH)), such that the 
Universe's homogeneity, isotropy, flatness, age and size problems are solved in the usual way j^]. 
Our inflationary model neither suffers from the usual runaway problem, presumed to be present in 
inflationary models with the Coleman- Weinberg potentials [3^, nor it suffers from the naturalness 
problem related to the usual fine tunning of the zero potential energy at the end of inflation. 

After the end of inflation, the field decays, the Universe reheats, and the potential energy is 
approximately given by and (jH^ . where H is now dominated by the radiation contribution, 
and can be viewed as a thermal (renormalisation) scale. The nature of the potential changes again 
in the matter era at a redshift, z 0.3, when dark energy starts dominating, and the Universe 
enters a new slow roll regime, driven by the effective potential ()4fjj) . While our model does not 
really solve the coincidence problem (why are the densities of dark energy and matter so similar 
today?), we point out that the required fermion mass ()61|) . whose fluctuations yield the correct 
amount of dark energy, is of the GUT scale, which is well motivated by particle physics. If during 
some earlier epoch in radiation or matter era the vacuum energy starts dominating (this may happen 
for example after a Higgs-like phase transition, or at the time of chiral condensate formation), then 
again becomes the relevant potential, driving dynamically Aefr toward smaller values. This means 
that phase transitions mediated by scalar fields or some other condensates, which may change the 
vacuum energy and thus the effective cosmological term, do not anymore pose a problem, since any 
vacuum energy gets eventually cornpensated. Even though it is not necessary to change the standard 
electroweak BEH mechanism j3] |2^, the Higgs potential may still happen to have a nonvanishing 
negative mass-squared term and a positive quartic self-coupling. A positive quartic coupling is not 
any more required, since the potential ultimately gets stabilised through the Yukawa couplings to the 
Standard Model fermions by the mechanism presented here. Since in this modified mass generation 
mechanism the Higgs mass would be of the order or smaller than the Hubble parameter today, 
tuh ~ -ffo ~ 10"'^^ eV, the Higgs particle would not be seen at the LHC or any future accelerator 
experiments. If the quartic self-coupling of the Higgs field is chosen to be strictly zero, then there 
is no hierarchy problem associated with large Higgs field radiative corrections (gauge fields radiative 



corrections still contribute). Furthermore, the issue of the electroweak vacuum stability |33|, and more 
generally the question whether this modified electroweak mass generation mechanism is consistent 
with all accelerator experiments, requires further investigation. 
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Finally, there are many possible improvements to this work: one can relax the assumption that 
the background space-time is (anti-)de Sitter, and work with quasi- ( ant i-)de Sitter space and more 
general expanding space-times; one can study how gradient corrections affect the results, etc. 
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